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TWO-BODY PROBLEM ON A SPHERE. REDUCTION, 
STOCHASTICITY, PERIODIC ORBITS 



We consider the problem of two interacting particles on a sphere. The potential of the interaction depends on the 
distance between the particles. The case of Newtonian-type potentials is studied in most detail. We reduce this system 
to a system with two degrees of freedom and give a number of remarkable periodic orbits. We also discuss integrability 
and stochastization of the motion. 



1. Introduction and historical notes 

The origins of the analysis of motion of mass particles and rigid bodies in constant curvature spaces 
(plane, two- and three-dimensional spheres, and Lobachevsky space) trace back to the classical XlXth 
century papers by Serret, Killing, Lipschitz, Liebmann, Schering, etc. One of the most comprehensive 
is the paper by W. Killing The author studied the motion of a particle in the "central field", a 
curved-space analog of the Newtonian field (the Kepler problem) . He also presented the generalization 
of the Euler two-center problem and iV-dimensional analogs of the above-mentioned problems. He 
derived and partially studied the equations of motion of a rigid body in these spaces and discussed 
some issues of the theory of the Newtonian potential. More than a century after that, the study 
of dynamics in constant curvature spaces attracted the interest of scientists again. This resulted in 
appearance of a series of papers, reviewed in [2]. Note, though, that some of the classical results were 
independently rediscovered in these papers. 

One of the major distinctive features of curved spaces is the absence of translational (Galilean) 
invariance, resulting in the absence of the center-of-inertia integrals and non-existence of the corre- 
sponding barycentric frame of reference. As a result, this takes us out of the realm of the classical 
celestial mechanics. 

For example, the two-body problem, when considered on a two-dimensional (or three-dimensional) 
sphere S 2 (S 3 ) or on a Lobachevsky plane (or space) L 2 (L 3 ), can no longer be reduced to the 
corresponding problem of motion in a central potential field (an analog of the Kepler problem) and, 
generally, is not integrable, as will be shown below. Hereinafter, we will discuss only the case of 
a two-dimensional sphere S 2 , though the obtained results can be easily extended to L 2 and their 
many-dimensional generalizations. 

First of all, let us note that the Kepler problem of motion in a central potential field on a sphere 
is integrable, and the orbits are ellipses with one of the foci at the center. An analog of the Newtonian 
potential on a sphere was first offered, apparently, by Serret in (1860), the same for a Lobachevsky 
plane was done by Lobachevsky himself and Bolyai. A generalization of the Bertrand theorem to a 
curved space was performed by Liebmann (1903). Analogs of Kepler's laws for S n and L n can be found 
in a number of papers, of which the most fundamental is the above-cited paper by W. Killing. Among 
the modern explorations of the Kepler problem from various perspectives, we recommend [HI El HI El- 
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So, an analog of the Newtonian potential on S 2 (L 2 ) is 

U = -~fctg6 (U = -jcth9), (1.1) 

where 9 is the longitude measured from the pole, at which the gravitating center is located, and 7 
is the gravitational constant. The potential (|l.lj) can be obtained as a centrosymmetric solution 
to the Laplace -Beltrami equations for 5 3 (L 3 ) or obtained from an analog of Bertrand theorem 
for S 2 (L 2 ) 0. 

In [2j, we state the restricted problem of two bodies on a sphere (or on L 2 ), where one of the 
bodies circumscribes a great circle on the sphere with constant velocity (i. e. moves freely along a 
geodesic), while the other moves in its field (with the potential (|1.1|) ) and does not affect the first 
body. In [S] (as well as in [2]), we offer also a numerical analysis of the Poincare map for this problem, 
which is not integrable, judging by its stochastic orbits; we also use the averaging method to study the 
perigee motion caused by the curvature of the space. In it is proved, for the restricted problem, 
the non-existence of an additional meromorphic analytic integral. 

2. Reduction of the two-body problem on S 2 

The general (not restricted) two-body problem on S 2 (L 2 ) was studied in respect to reduction in [T2] . 
where the quantization problems were also discussed. The general method of reduction, offered in |12j . 
was first suggested by E. Cartan and further developed by J. Marsden and A. Weinstein. This method 
is based on rather formal argumentation, which, at the same time, has a profound differential geometric 
interpretation. Unfortunately, the method does not always gives and suitable expressions, which have 
been sought after in the classical celestial mechanics (for example, when reducing the three-body 
problem). 

In this paper, to reduce the unrestricted two-body problem on a sphere, we use Bour's reduc- 
tion [Tn| of the spatial three-body problem (similarly, Radau's reduction can be used JB]). In a sense, 
the reduction we need can be regarded as its special case. Connection between the mentioned problems 
is based on a famous observation by Jacobi: once a barycentric frame of reference was introduced, 
which is equivalent to reduction to the center of inertia, the classical three-body problem could be 
reduced to the problem of two mass points moving in the fixed space Oxyz with a potential, which 
was a function of the points' distances from the fixed center O and of the angle between their radius- 
vectors. Such a reduced system has a (vector) integral of angular momentum, and Bour's reduction is 
an effective reduction of the system's order using this integral (elimination of a node). With that, two 
more degrees of freedom are eliminated. In the case of a sphere, the distances r±, r% between two points 
and the center are fixed, so in the Bour-reduced system we only have to put r\ = r<i = R, f\ = r<i = 
{R is the sphere's radius). The specified reduction is valid for the potential, which arbitrarily depends 
on the distance between the points (the two-body problem is sometimes referenced as E. I. Kugushev's 
problem pQ). 

Below we give the calculations in more detail. Let Oxyz be a fixed frame of reference, and O^rjQ - 
a moving frame of reference, such that the plane Or/^ contains the mass points mi and m-i- The axis rj, 
being an intersection of the plane (77^) with (xy), is the line of nodes (see Fig. The mass points 
are described by the coordinates 6\ and 62 (the angular polar coordinates in the plane Or]^) and (p, ip 
(these angles specify the location of the frame 0£r/£ relative to Oxyz). It is easy to see that all these 
variables are identical with the Euler angles from the rigid body dynamics (this analogy was first 
noted by J. Silvestr). 

In terms of the introduced variables, the Cartesian coordinates read: 

Xi = —R(cos <p cos ip sin 6i + sin tp cos 9i), 

yi = —R(cos(psmipsin9i — cosipcos9i), i = l,2. (2-1) 
Zi = R sin 9i sin tp, 
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Fig. 1. 

Using H2.1j) . the Lagrangian function of the system can be written as 

Jz? = y(mi(fli +-0cos<^) 2 + m 2 {9 2 + ^cosp) 2 + Jp 2 + 

+ R sin 2 pip 2 - 2Lpip sin ip) - U{Qx - 2 j) , (2.2) 

where 

J = mi sin 2 9\ + m 2 sin 2 6 2 , K = mi cos 2 8\ + m 2 cos 2 02 > 
L = mi sin #i cos 6\ + m2 sin 6 2 cos #2- 

Let the axis Oz be directed along the vector of angular momentum \M\ = c. The Lagrangian 
equations with the Lagrangian ()2.2j) admit the integral 



- M z = c = const (2.3) 



dip 

and invariant relations 

ff = M " = °' ff = M -^ = °- (2 " 4) 

Indeed, 

M„ = M zxv ijj, M zxr) = -M v ijj, 

where M z , M„, M zxr) are the projections of the angular momentum M onto the given axes ^H] (see 
Fig. 1). The invariant manifold ()2.3|) . (|2.4|l is explicitly given by: 

• c j . cL sin 

i J = ^ . o,„ rr» ^~ 



mim 2 R 2 sin 2 (6>i - 6> 2 ) mim 2 R 2 sin 2 (6>i - 9 2 ) 

(2.5) 

mi 0i + m 2 6>2 

COS ip = r. 

c/R 2 — (mi + m 2 )ip 

It follows from (|2.3fl that the variable ^ is cyclic, so we can reduce the system's order by one 
degree of freedom. Besides, the variable p is also cyclic on the invariant manifold (|2.4[) . and here we 
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again can reduce the number of degrees of freedom by one, using the Routh reduction procedure. The 
Routh function is 



St = J? - - ^ = ^-( mi 0l + m 2 9 2 2 ) - 

dtp 9 <p 2 

R 2 J(m 1 e\+m 2 e 2 ) 2 c 2j 



(2.6) 



2(J(mi+m 2 )-mim 2 sin 2 (6>i-#2)) 2mim 2 -R 2 sin 2 (6q-6>2) 
Using the Legendre transformations 

d& r>2a rriiR 2 J {m^ + m 2 2 ) . 

Pi = — — = rriiR Vi ^ , i = 1, 2, (2.7) 

<9#j J(mi + m 2 ) — mim2 shr (#i - 6*2) 

we obtain the equations of motion of the Hamiltonian form with the Hamiltonian 

2 2 

* = (w + m + (c2 " (pi + ^ 2)2) ^ — zt) + u ^ ~ (2.8) 

2R Z V 1 12 mim 2 sm z (6'i - 6> 2 ) J 

This is the sought-for reduced system with two degrees of freedom. 

The absolute motion (i. e. motion in the fixed frame of reference) can be obtained, using the two 
additional quadratures given by (|2.5|) . 

Remark. On the manifold (|2.5|) , the domain of possible values of the new variables pi and pi is bounded 
because it follows from 12.5fl that cos ip = ^ 1 -7 ^ 2 and | cos^l ^ 1. 

So, the phase space of the reduced system is parametrized with the canonical vari- 
ables 9 1 , 6 2 , Pi, p 2 - The angle variables 1 , 6 2 mod 2tt define a two-dimensional torus (configura- 
tion space of the reduced system), while the conjugate momenta p 1 , p 2 form the strip \p l + p 2 \ = 
= |ccosy| ^ |c|. Besides, it is necessary to identify the points of the phase space with coordinates 

(#1, #2) Pi, P2) and ^1' Pi) P2Y ( 2 - 9 ) 

Indeed, according to (|2.1|) . the different sets of coordinates, (0^ 2 > V 3 ) an d (vr — 9 1 ,ir — 8 2 , it — 

— <Pi, Tr + tp), define the same configuration of the bodies on the sphere. For the velocities we have (ir — 

— = (tt — 2 )' = — # 2 an d 5 using (|2.7|) . obtain (|2.9|) . We believe it would be an interesting to 
develop a better geometric description of the reduced phase space. 

The Hamiltonian Q2.8[) is a homogeneous, but not positively defined, quadratic function of the 
momenta. Since the system (|2.8j) is not natural, well-developed methods of topological and qualitative 
analysis j5] can hardly be applied to it. 



3. Particular solutions and the Smale diagram of the two-body problem on S 2 

Let us start with the general problem of n interacting bodies. The potential of interaction U depends 
on the distances between the bodies. 

A particular solution to the n-body problem will be called a rigid-body motion if the distances 
between the bodies remain constant. A particular solution to the n-body problem will be called a 
stationary configuration (relative equilibrium) if all the bodies rotate uniformly, with the same angular 
velocity, about the axis containing the vector M. 

It is obvious that any stationary configuration is a rigid-body motion, but the reverse is, generally, 
not true. 
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3.1. Stationary configurations and the Smale diagram 

The stationary configurations in the general n-body problem on S 2 can be obtained as the critical 
points of the reduced potential. 

Proposition 1. The stationary configurations (relative equilibriuma) of the n-body problem cor- 
respond to the critical points of the reduced (effective) potential 

C/, = [/-iu, 2 /=[/-^, (3.1) 

where u is the angular velocity of the configuration, / is the total moment of inertia of the system 
with respect to the axis of rotation, M 2 is the squared total angular momentum of the system of 
bodies in a fixed frame of reference. 

Proof. 

Consider in the frame of reference, rotating with the angular velocity oo about the axis z. In this 
frame, the Lagrangian (in terms of spherical coordinates) reads: 

Se = \ m^ie 2 + sin 2 9^\) + u ]T m^R 2 sin 2 9, - U + ± W 2 /, (3.2) 

where I = ^2 m i R 2 sin 2 9 i is the total moment of inertia with respect to the z-axis (i is the number of 
a particle) . 

Then, using equations of motion with the Lagrangian (|3,2j) and bearing in mind that 9 i = <p i = 
(stationary solution), we obtain the following conditions of the relative equilibrium: 

4( c/ "^ 2/ ) =0 ' i=1 '-' w - (3 - 3) 

To finish the proof, we have yet to express the angular momentum of the stationary configuration 
in the fixed frame of reference: 

M x = M y = 0, M Z = Y, m.R 2 sin 2 Ofa =^ ^ ^ *i = u 

Now we consider the case of the generalized Newtonian potential and give a complete description 
of the relative equilibriuma of the two-body problem on S 2 . According to Proposition^ it is necessary 
to find the critical points of the reduced potential 

2 2 2/2 2 \ 

U* = —asR ctg 9 12 — 2^ R \ m i s ^ n @i + m 2 s i n &2J > 

where 9\2 is the angular distance the particles on the sphere, and as is the gravitational constant. 

It can be shown that for relative equilibriuma of two particles on S 2 , the following two observations 
are true: 

1. While moving, the particles never leave the great circle that contains the axis of rotation. 

2. When a? > (the case of attraction) the particles stay on the great circle, with the axis of 
rotation between them (see Fig. |2J). 

Indeed, the cosine of the angle between the mass points can be rewritten as 

£ = cos 9 12 = cos 9 l cos 6*2 + sin 6^ sin #2 cos(<^i — ^2), 

and we find that 

8U* dU* seR 2 m im2 s n 

— — = — — — = — sin fri sin 9 sm Wi — gj 9 = U. 

di Pl dif 2 (1- £2)3/2 1 2 Wl Vl) 
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dU 

Hence, either tp-, = (f2, or ip 2 = ip-i + ix. If ipi = (f2 the equations = do not have a common 
solution. 

So, to find stationary configurations, it is necessary to find the critical points of the function 
J7* = —seR 2 m 1 m 2 ctg{6 l + 9 2 ) — ^oJ 2 R 2 (m 1 sin 2 9 1 + m 2 sin 2 9 2 ) 
where the angles 6 1 and 6 2 are measured from the axis of rotation (Fig. |2J) . These angles are given by 

2 aSTOg 2 33777^ 

ijj sin 0i cos 9\ = 0, u sin # 9 cos 6*2 = 0. (3.4) 

1 sin 2 (^ + 02) sin 2 (0 1 + 2 ) V ; 



n 


Mi 

>./ 


\?7t| 
A// \ 













Fig. 2. 

From this, we obtain the equation that generalizes the planar "law of the lever" and relates the 
angles 9 X and 9 2 

m l sin 29 l = m 2 sin 29 2 . (3-5) 

Using (|3.4j) and ()3.5j) . one can show that for 33 > 0, either 6\, 9 2 < tt/2, or 9±, 9 2 > tt/2. These 
configurations are identical up to the reflection about the equatorial plane. Therefore, we will assume 
that 

< 9 X , 9 2 < f . (3.6) 
If we put m 1 < 77i 2 , then the equation has two roots on the interval (|3.6|) 

fl (l) 1 ■ f m ± ■ on\ fl( 2 ) 7T 1 ■ f m l ■ nn\ (o 7 N 

^2 =2 arcsm I sm 28 1 1 , 2 2 — 2 arcsm \ 1 ) ' ' ) 

here 0^ < X , and 9 { 2 ] > 9 X . 

The first solution 9 2 describes a configuration, where the heavier particle moves along the smaller 
circle. As R — ► 00, this solution tends to ordinary stationary configuration of the two-body problem 
on a plane. 

(2) 

The second solution 9 2 describes a configuration, where the heavier particle moves along the 
larger circle; this solution does not have an analog in the two-body problem on a plane. 

Let us plot a bifurcation diagram (see Jl]) for the two-body problem on S 2 on the integral 
plane (c 2 = M 2 /R 4 , E = h/R 2 ), where h is the energy integral (Fig. |3J). When m 1 7^ m 2 , we obtain 
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two curves, corresponding to the solutions of Q3.7JI (see Fig. When m 1 = m 2 , the curves merge 
at a point P (Fig. |3t>) ; in this case, also according to (|3.7|) . we have either 9 t = 9 2 , or 9 1 + 2 = n/2. 

Here we also give a bifurcation diagram for the two-body problem on a sphere in the presence 
of a potential that, unlike the gravitational potential, has no singularities (see Fig. 0J). For that 
potential, in the paper [Q a geometrical (topological) analysis of the integral manifolds M 2 = const, 
h = const and of the level surface of the reduced system's energy integral is performed. Note that, 
due to singularities at the poles in the case of the gravitational potential it is impossible to 

directly apply the results of [Q to this case. So, the complete topological analysis of three-dimensional 
isoenergetic manifolds does not seem to have been done. 




Fig. 3. Bifurcation diagram for the two-body problem with the gravitational potential U = —sem 1 m 2 ctg# 




Fig. 4. Bifurcation diagram for the two-body problem on a sphere with the potential U = aem 1 m 2 (l — cos 9) 



3.2. A particular solution in the case of equal masses 

The ra-body problem on a plane has a particular solution (collinear configuration), in which all the 
particles move along ellipses, while staying at all times on the same straight line. This solution was first 
offered by Euler (for the three-body problem) and Moulton (for the n-body case) [5]. In [5], Moulton 
also proved that for the case when all masses are different the number of such collinear configurations 
(for different permutations of the particles) is equal to Let us show that analogous (pulsating) 
solutions on a sphere do not, generally, exist. 

Proposition 2. For the two-body (n-body) problem on S 2 , when m i ^ rrij, there are no solutions 
other than a stationary configuration, in which all the particles stay on the same great circle that 
contains the axis of rotation. 
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Proof. 

Let Oz be the axis of rotation. The equations of motion in spherical coordinates are 

d 8U ( 3-8 ) 

— (m 8 sin 2 = -— = V UL sin9 i sin^ sinfo - p.), 

where ^ • = cos Q i cos ■ + cos 9 i sin ■ cos(<^ — <Pj)- 

If the particles always stay on a great circles that contains the axis of rotation, then either ip i = <p-, 
or <fj = ip i + 7r and, besides, 

<Pi = l Pj = ' ) P, h j = 1) • • • , 
Therefore, according to 1)3.8)1 . the following must hold: 

m i sin 2 0^ = m,jC 2 , q = const (3-9) 

Hence, all the angles 0j depend only on the constants q and the function A 2 (t) = l/tp(t) 

C i = sm6 i = c i \(t). (3.10) 

Differentiating this and substituting into the first equation of ()3.8|) . we obtain 

?7 V» (< , c 2 AA 2 \ 3 r — qu 

, A H — - = rrijCjA a d 1 - c z A 2 - (3.11) 

^1 - Ci A 2 \ i-c 2 A 2 y ' v 4 ^ 

where tttt = t^t-^tt = a/1 — c 2 A 2 -^— . The required solution exists if (13.111) can be made independent 

36, dQde, V dd 

on the particle's number i. Provided that A ^ 0, this can be done if , d i are assumed to be the critical 
points of the reduced potential (|3.1|) and all c { are assumed to be equal, i. e. 

= uj 2 m i sm6 i cos9 i , q = c = const. (3.12) 

Then, taking into account (|3.1U|) and 1)3.11(1 . incompatible systems are obtained 

mj sin 29 i = ±nij sin 29 j and sin 9^ = ± sin ■ , 

which does not have a solution m- ^ m-. ■ 

Let us consider more closely the case of equal masses in the two-body problem, m 1 = m 2 = 1 
(with ip 2 = <fi + vr and 9 1 = 2 = 0). Using directly 1)3.8(1 and ()3.9() . we obtain 

1 p = ( p 1 = i P = c = sin0cos0V> 2 rfrr, c = const. (3.13) 

sir sir 29 v 7 

These equations describe the motion of a particle on a sphere in the field of an attracting center with 
the potential 

U = -Iaectg20. 

For this potential, the orbits are unclosed, which showes another difference between curved and 
flat spaces. It is known that, due to uniformity condition, the particles' orbits in flat spaces are ellipses 
with one of the foci at the center of mass. 

Nevertheless, orbits on a sphere are closed in some rotating frame of reference. Moreover, one can 
choose the velocity of rotation so that both particles will move along one and the same curve (such 
motion is called a relative choreography see Fig. O B±, B2, B3). 
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3.3. Rigid-body motions in the two-body problem 

Consider the motion of particles in a flat space M 3 in a fixed frame of reference. In this case, the relative 
equilibriuma of the two-body (n-body) problem in the center-of-mass reference frame, generally (i. e. 
when the velocity of the center of mass is not zero), correspond to the simplest rigid-body solutions. 
The mutual distances between the particles are always the same. We show that in a curved space 
(on if? 2 ), lacking the notion of the center of mass, all the rigid-body solutions of the two-body problem 
are stationary configurations discussed above (|3.4|) . Q3.7JI . Therefore, the corresponding family of 
configurations has less free parameters than such a family in a flat space. 

Theorem. When the interaction potential depends only on the mutual distances, all rigid-body 
solutions of the two-body problem on S 2 , are stationary configurations. 

Proof. 

The proof is based on successive differentiation of the invariant relation 

(x, y) = cos 6*o = const, (3-14) 

where x, y are the radius- vectors of the particles on S 2 in R 3 , and finding the maximal invariant 
manifold that contains these motions. 

Now let us perform an algebraic reduction of our system by restricting it onto a level surface 
of M. For the variables 

M = vtl x x x x + rn 2 y x y, L = fi(x x x — y x y), 

wfiere a = ; , we fiave nine equations: 

^ m 1 + m 2 

L = U'^x x y 

x = x x M + -^-x x L, (oik\ 

m 1 + m 2 rn-i (3.15) 



1 1 
y = i y x M - —-y x L. 



Here U(£) is the potential energy of the interaction, expressed in terms of the cosine of the angle 
between the particles £ = (x, y), and vector M is a parameter. Due to the condition (x, x) = 
= {Vi y) = and definition of x, y, the following relations hold: 

x =y =1; 

'x, M) + {x, L) = 0, -^ 2 -(y,M)-(y,L) = 0. ^ 



m 1 + m 2 m 1 + m 2 

The energy integral in this case reads 



E = 4-L 2 + U(0 = const. 



Differentiating the relation (x, y) = const by time and using 1)3.15(1 and (|3.16|) . we get 

(x, y x L) = 0, 

±L(x, y) + U'(£){x x yf + g(x, M)(y, M) = 0, (3.17) 
fx, M x L) = 0. 
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It follows from the second and the fourth equations of 1)3.17)1 that the vectors x, y, M, and L 
belong to the same plane. Express M, L in terms of x, y: 

% , (M, x)- cos 6 (M, y) (M, y) — cos 9 (M, x) 

M = ^T~a X + ^ y ' 

sin f/ n sin a n 

(3.18) 

(L, x) - cos 9 (L, y) (L, y) - cos9 (L, x) 

L= ~z7> X + ^T2 y; 

sm 9 sm 9 

substituting this into 1)3.15)1 and 1)3.17)1 . we have (M , x) = const, and (M, y) = const. Thus, the 
particles move in circles about some axis defined by the vector M, and during the motion the particles 
and the axis always lie in the same plane. ■ 

3.4. The case of M = and the two-body problem on a circle. Collisions 

In the two-body problem on a sphere (with an arbitrary potential, depending only on the inter-particle 
distance), the motion with zero total momentum {M = 0) is a special case. We will show that the 
particles in this case move along a great circle, which holds its position in a fixed frame of reference. 
Indeed, for M = 0, equations ()3.15)) read 

± = lk[ xxL > y = ~Wqy xL > L=U'^xxy, (3.19) 
while the relations ()3.16j) take the form 

(x, L) = (y, L) = 0. (3.20) 

Calculating the orbital derivative of the normal to the particles' plane n = X X — under the 
\x x y\ 

flow of 1)3.19)1 and using l)3.2U)) . we get h = 0. Thus, the plane that contains the radius- vectors of the 

particles is fixed, and the particles move along a great circle. 

Remark. In the two-body problem on a plane, when the momentum is zero, the particles move along the 
same straight line, fixed in the center-of-mass reference frame. This frame moves uniformly and in a straight 
line. 

Since for zero momentum M = 0, the two-body problem on a sphere is reduced to a system on 
a circle, we will consider the latter problem in more detail (for arbitrary M). We will also show that 
in this case there is a center-of-mass reference frame, in which the behavior of the particles does not 
depend on the value of M. The new variables are 

V>= 1 1 , — , = 0j -0 2 , 3.21 

where X and 2 are the angles that define the positions of the particles on the circle. So, the Lagrangian 
reads 



££= U mi +m 2 )i; 2 + ^9 2 -U(9), ti 



2 V 1 2^ K n r m l + m 2 ' 

This implies that the "center-of-mass angle" ip varies uniformly with time: 



; \ M \ , 

^ = l^T 7 = const - ( 3 - 22 ) 

H {m 1 + m 2 ) 



The time dependency of is given in terms of quadratures 

d6 = 

^-\h l -U{9)) ~ 

where O , t , and ^ are some constants. 



t-t , (3.23) 
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For the generalized Newtonian potential U(9) = —sem 1 m 2 ctg9, we find that the particles will 
collide in finite time. If t = at the time of collision, the time dependency of the angle can be written 
in the form of the Puiseux series 

oo 

6(t) =t 2 / 3 J2 c J n/3 - (3-24) 

n=0 

Here c 7^ 0, while the coefficients c n with odd indices are zero (c2 n +i =0). In other words, 9(t) is 
an even function of time. Consequently, the particles undergo perfectly elastic collision (for a suitable 
regularization). 

Remark 1. The proof of (|3.24|l is based on the analysis of (|3.23() near 8 = 0. Introducing a new variable x = 
= tg0, we express <|3.23[1 in the following form (for t =0): 

s/xdx I " ~ h 1 



C, C = \j 238(77^ + m 2 )t, a 



(1 +x 2 )Vl + ax 

U 

Using the Taylor expansion of the integrand, we obtain the following series equations defining x(t): 

00 



n=0 

where g n are constants that depend on a. After eliminating u, we have 



*(0 = C 2/3 £ ^ 



/3 

'nS 

n=0 



Remark 2. We believe it would be interesting to generalize the classical results obtained by Zigel, concerning 
zero measure of collision orbits, to this problem (i. e. to specify the collision manifold). 

Let us also mention a long-discussed question [21 [7] about a possibility of extending Sundman's results (i. e. 
solutions in power series) to the two- and three-body problems in curved spaces. Unfortunately, the answer to 
this question has not yet been found and is, most likely, negative. 



4. The Poincare section. Numerical analysis 
4.1. The surface of section. Chaos 

To study the system's behavior numerically, we plot the Poincare section of the reduced system (|2.8jl in 
the following way. The two-dimensional surface of section of the three-dimensional energy level Jif = 
= E is given by p 2 = const. Unlike one-and-a-half degree of freedom systems (or trivial two-degree 
systems), the surfaces of section in our case might have a complicated structure. Moreover, the 
projection of such a surface on a plane (for example, (p t , 9^) has singularities. To avoid them, we will 
plot these surfaces in the three-dimensional space of variables (6 1 , p v 9 2 )- Some examples of surfaces 
of section p 2 = are given in Fig. [21 The section's type depends on the location of the point (E, c) on 
the bifurcation diagram (Fig. EJ). The phase flow of the system (|2.8|) generates a point Poincare map 
on these surfaces. The Poincare maps for the case of equal masses (for some values of c and E) are 
shown in Fig. 03 According to 1)2 .9|) the points with coordinates (9 ± , 9 2 , Pi) and (it — 9 X , tt — 9 2 , ~Pi) 
should, generally, be identified. For large energies, there are chaotic motions on the sections (Fig. 
f, 1), which prevent the existence of an additional analytic integral of the system (|2.8|) . 

Remark. Non-existence of an additional analytic integral for (j2.8(l can also be proved using Poincare's 
methods |10j . i. e. analyzing the secular terms of the perturbed system and the corresponding number of non- 
degenerate periodic orbits. In this connection, the expansions of the disturbing function, given in |1U| . may be 
useful. These expansions were used in the proof of non-integrability of the classical three-body problem. 
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Fig. 5. Energy level sections (a, c, e, g, k) and corresponding Poincare maps (b, d, f. h. 1) for different values 
of E and c. 

4.2. Periodic solutions and choreographies 

An important part in the classical celestial mechanics is played by periodic orbits in a center-of-mass 
reference frame. Such orbits were studied in many papers [15] . Recently, a new class of periodic 
solutions in the problem of n equal-mass bodies has been discovered. These are so-called (relative) 
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Fig. 5. (Continuation.) Energy level sections (a, c, e, g, k) and corresponding Poincarc maps (b. d, f, h, 1) for 
different values of E and c. 

choreographies — solutions, when the particles follow each other (in a rotating frame of reference) 
along the same curve, with time shift T/n (^Bj 5 see also [3]). 

A similar part in the two-body problem on a sphere is played by periodic solutions of the reduced 
system (|2.8|) . The simplest periodic solution that can be expressed analytically is (|3,13|) . There are 
other similar solutions, which, however, cannot be presented by analytical formulas. A suitable tool 
to analyze these solutions is the Poincare map constructed above. Fixed points of the least order 
correspond to the simplest periodic solutions (see Fig. EJ). These solutions define choreographies of 
particles (the particles move along the same curve) in some rotating frame of reference [3J E3 • F° r 
the least-period periodic solutions, the projections of such choreographies on a plane, perpendicular 
to the axis of rotation, are given in Fig. H3 Figures A±-Aq and B\-B^ show two different stable 
solutions, while figures C\—Ci show the unstable solution of the reduced system for different energies. 
Figures B 1 , B 2 , and B 3 show the analytic solution (|3.13|) . 

Remark. It can be shown, using the methods of [3j , that for a fixed c, there exists a countable set of values 
of E, for which the particles form an absolute choreography, i. e. move along a closed curve in a fixed space. 

Remark. When m 1 ^ m 2 , choreographies are broken — each particle moves along its own closed curve (in 
a rotating frame of reference). 
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Fig. 6. Relative choreographies for the problem of two equal-mass bodies (mi = ni2 = 1) on a sphere, c = 1, 
fl= 1. 
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